arXiv:1508.01699v3 [physics.optics] 15 Apr 2016 


Point Interactions, Metamaterials, and PT-Symmetry 


Ali Mostafazadeh* 

Departments of Mathematics and Physics, Kog University, 
34450 Sariyer, Istanbul, Turkey 


Abstract 

We express the boundary conditions for TE and TM waves at the interfaces of an infinite 
planar slab of homogeneous metamaterial as certain point interactions and use them to com¬ 
pute the transfer matrix of the system. This allows us to demonstrate the omnidirectional 
reflectionlessness of Veselago’s slab for waves of arbitrary wavelength, reveal the translational 
and reflection symmetry of this slab, explore the laser threshold condition and coherent per¬ 
fect absorption for active negative-index metamaterials, introduce a point interaction modeling 
phase-conjugation, determine the corresponding antilinear transfer matrix, and offer a simple 
proof of the equivalence of Veselago’s slab with a pair of parallel phase-conjugating plates. We 
also study the connection between certain optical setups involving metamaterials and a class 
of PT-symmetric quantum systems defined on wedge-shape contours in the complex plane. 
This provides a physical interpretation for the latter. 

PACS numbers: 03.65.Nk, 78.67.Pt, 42.25.Bs, 42.65.Hw 

Keywords: Metamaterial slab, point interactions, invisibility, laser threshold condition, phase- 
conjugation, PT-Symmetry. 


1 Introduction 

Material with negative permittivity and permeability have nnnsnal electromagnetic properties [1]. 
This observation was originally made by Veselago in 1968 [2] , bnt it was not nntil the year 2000 that 
it began to receive the attention of other physicists and become a focus of intensive research EH]. 
The interest in this type of metamaterials was triggered and reenforced by two major developments. 
First, they were no long hypothetical substances but could be realized in laboratories [5]. Second, 
they were shown to have remarkable applications in producing perfect lenses [3] and invisibility 
cloaks |1]. 

The purpose of the present article is to explore some of the basic properties of metamaterials 
using the machinery of point interactions and transfer matrices. Specihcally, we consider the scat¬ 
tering problem for transverse electric (TE) and transverse magnetic (TM) waves interacting with 
an inhnite planar slab of homogeneous metamaterial. We identify the effect of the electromagnetic 
boundary conditions with the presence of certain point interactions and use them to obtain an 
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Figure 1: (Color online) A wedge-shaped contour in complex plane. 


explicit expression for the transfer matrix of the slab. This in turn allows us to investigate the 
invisibility properties and symmetries of the Veselago slab and derive explicit expressions for the 
laser threshold PE] and antilasing conditions ili for active metamaterial slabs. Similarly we 
address the problem of modeling the phenomenon of phase conjugation [ 10 ] using certain antilinear 
point interactions and provide a very simple description of the equivalence of a Veselago slab to a 
pair of parallel phase-conjugating plates m- Next, we examine a recently studied spectral problem 
[TS] that arises as a simplification of a metamaterial cloaking ring model [13]. We offer a physical 
interpretation of this problem in terms of the transverse modes of a rectangular waveguide half- 
hlled with a metamaterial with permittivity —1 and permeability -|-1, and establish its relation to 
the spectral problem for a PT-symmetric inhnite square well potential defined on a wedge-shaped 
contour A of the type shown Fig. [TJ We can express it in the form 

r 0 for |Re(j)| < «. 

( cx) otherwise, 

where 3 marks a point on A and £ is a positive real parameter. The quantum system associated 
with this potential was studied a decade ago in the context of obtaining a real description of VT- 
symmetric systems defined on complex contours [H]. This was a crucial step towards the application 
of the formalism of pseudo-Hermitian quantum mechanics [13] for such systems. It was also the 
key to the subsequent results on the relation between PT-symmetric wrong-sign quartic potential, 
v{x) = and the usual quartic anharmonic potential, v{x) = x^ + ax"^, [16]. Here we relate 
PP-symmetric inhnite square well potentials dehned on wedge-shaped contours to the physics of 
metamaterials. 


2 Metamaterial Slabs and Point Interactions 

Consider an inhnite planar slab S of thickness L that is hlled with a homogeneous (meta)material 
and placed in vacuum. If we take a cartesian coordinate system where the boundaries of S corre¬ 
sponds to the z = 0 and z = L planes, we can write the permittivity e and permeability fi of this 
system according to 

_ f eeo for ze[0,L], _f ppo for ze[0,L\, 

^ \ £0 for z^[0,L], ^ \ Ho for z^[0,L], 

where eo and e are the permittivity of the vacuum and the relative permittivity of the slab, and ho 
and H are the permeability of the vacuum and the relative permittivity of the slab, respectively. 
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Figure 2: (Color online) Graphs of the TE (on the left) and TM (on the right) waves incident on a 
slab of (meta)material (the colored region). 


Now, consider a time-harmonic TE (respectively TM) wave incident on S such that the electric 
field E of the TE (magnetic field H of the TM) wave lies along the y-ax.\s and the wavevector k 
makes an angle 9 with the 2 ;-axis as shown in Fig. |2l Then, 

k = kxCx + bzCz, kx'-=ksiYi6, kz'.= k cos 9, (3) 

where k := \k\ = oj/c is the wavenumber, oj is the angular frequency of the wave, c := 1 / ^Eq po is 
the speed of light in vacuum, and e^,, Cy, and are the unit vectors along the positive x-, y-, and 
z-axes, respectively. 

With these conventions in place, we can write the TE and TM solutions of the Maxwell equations, 

V-(£E) = 0, V-(/i#) = 0, (4) 

edtE = VxH, ydtH = -V X E, (5) 


as 


E = <^{z)ey, 

H := [ik-^^'{z)ex + sm9^{z)ez] , 

E := [ik-^jr'{z)ex + sin 9 J^{z)ez] , 

H = 

Here S’ and JS are solutions of the Helmholtz equation. 


TE ; 


TM : 


-^"(z) -I- k'^{n^ — sin^ 9)i/j{z) = 0, z ^ {0, L}, 


n := e fic , 


which satisfy the matching conditions, 

Siz^)=Siz^), 

JS{Z+) = je{zo), 


S'iz^)=aJ^S\z^), 


( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 
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at zq = 0, L, for every function / : M —)■ M, stands for the limit of f{z) as z tends to zq 

from the left/right, and 

TE ) 

For the system we consider, 


a™ := 


g(4) 

e(Zn)' 


( e ft for z ^ [0, L], 
\ 1 for z ^ [0, L], 


On 


, TE 



^TM 

Oq 



( 11 ) 


The term metamaterial refers to the unusual situations where the real part of i, /t, or both have a 
negative sign. The best known example is a Veselago metamaterial, where e = /i = —1, [2]. As seen 
from the above equations, all the interesting properties of the latter follows from the corresponding 
matching conditions (E]) and fIlOp at the boundaries. This observation has motivated the authors of 
m to argue that the same effects can be realized using a pair of parallel phase-conjugating plates 
placed in vacuum. In the remainder of this section we consider a more general situation where i 
and fi are arbitrary complex numbers with a negative real part and try to construct an equivalent 
system consisting of a slab S of ordinary material and a pair of plates that realize the relevant 
boundary conditions in such a way that the permittivity e and permeability fi of S differ from those 
of <S only by the sign of the real part of e and fi, i.e.. 


Re(e) = —Re(£), 
Im(£) = Im(d), 


Re(/i) = —Re(/t), 
Im(/i) = Im(/t). 


( 12 ) 


Let us view the Helmholtz equation ([8]) as the time-independent Schrodinger equation 



- <(z) + v{z)'ip{z) = klij{z), 

(13) 

for the potential 


v{z) := k‘^{l — n^) = fc^(l — n^). 

(14) 

where 


h := ± sec — sin^ 9, 

(15) 


and the undetermined sign in this relation is fixed by the requirement that the real part of h and n 
have identical sign, i.e., 

Re(n)Re(n) > 0. (16) 


Then the presence of the above-noted plates corresponds to the addition of certain point interactions 
to the potential (ITTp . Notice also that (ITTp is, in general, a barrier potential that vanishes outside 
[0, L] and has a (possibly) complex hight fc^(l — ifi) . 

The determination of the point interactions that map the boundary conditions for the metama- 
trial slab S to those of S requires the knowledge of the effect of the addition of a point interaction 
to a hnite-range scattering potential. This is most conveniently acquired using the transfer matrix 
formulation of scattering theory. 
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3 Transfer Matrix for a Metamaterial Slab 


For a scattering potential v{z) which decays sufficiently rapidly as z ^ ±cxo, every solution of the 
Schrodinger equation (IT^ has the following asymptotic form. 

ij{x) for z ^ ±cx>. 


The transfer matrix of v{z) is by dehnition the 2x2 matrix M satisfying 


B+ 


M 


A. 


Let us recall that the left- and right-incident scattering solutions of flT^ are dehned by the asymp¬ 
totic boundary conditions: 


iJiiz) iV, X 1 

„ikzZ 1 jDl„—ikzZ 

for 

2; —)> —00, 

( 17 ) 


for 

z ^ 00, 

'lpr{z) Nr X < 

f 'J'r^-ikzZ 

for 

z —00, 

( 18 ) 

i Q-ikzZ _j_ J^r^ikzZ 

for 

z ^ 00, 


where Ni/^ are normalization constants, and and T*/'’ are respectively the left/right reflection 
and transmission amplitudes. These are related to the entries and determinant of M according to 




M21 

W2' 


RT 


Mi2 

W2' 


det M 
M22 ’ 



(19) 


Using the Wronskian identities satished by the solutions of flT^ . one can show that detM = 1, 
[mill]. This relation holds for both real and complex scattering potentials. It implies that the 
left and right transmission amplitudes are equal. Because of this we use T for T*/'’, whenever we 
consider a scattering problem obeying det M = 1. 

An important property of transfer matrices is their composition rule. This states that if v{z) is 
the sum of two potentials V-{z) and v+{z) such that the support of V- lies to the left of that of n+, 
i.e., there is some a G M, w_(z) = 0 for z > a and v^{z) = 0 for z < a, then the transfer matrix of 
v{z) has the form M = M+M_, where M± is the transfer matrix of v±{z). 

The above description of the transfer matrices, scattering solutions, and reflection and transmis¬ 
sion amplitudes applies also for the scattering due to a point interaction. The latter is a zero-range 
interaction whose effect is to impose matching conditions for 'ijj{z) and ip^z) at a discrete set of 
points Zq in R. For example consider a point interaction with a single center zq that is dehned by 


-1 

= B, 

1- 

1 0 

1_ 

1- 

■^ZQ 

1- 

° 1 


( 20 ) 


where is an arbitrary invertible matching matrix [/|. Then for z ^ zq the system corresponds to 
a free particle and we can easily determine its transfer matrix. This gives [19] 




^There are also other types of point interaction. See for example those considered in [20] . 


( 21 ) 
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where 


N.o:= 


^ikzZQ g ikzZQ 

n ]/^ ^^kz^Q rt c> ikz^Q 

crvzC' LrvrC' 


( 22 ) 


An immediate consequence of (12T1) is the existence of anomalous point interactions that violate the 
condition det M = 1. Clearly they correspond to matrices that do not have a unit determinant 


We can identify the electromagnetic boundary conditions (jH]) and fllOl) with a pair of point 
interactions with center zo = 0,L and diagonal matching matrices: 


gTE/TM 


1 0 

n TE/TM 
0 azo 


Substituting this relation in fl2T]) . we can express the corresponding transfer matrix as 


M™/tm _ 


(23) 


(24) 


where 

1 + a (1 - 
(1 - 1 + a 

Clearly, det = det Therefore, electromagnetic boundary conditions 

generally correspond to anomalous point interactions. However, often these come in pairs in such 
a way that their contribution to the determinant of the transfer matrix of the system as a whole 
cancel and we again hnd a total transfer matrix that has a unit determinant. 

Next, we conhne our attention to the TE waves. Let denote the transfer matrix for the 

TE waves and M* be its value for a hctitious slab iS* of the same size and position that is hlled 
with material with permittivity e* := eoifi and permeability fi-), := fiQ. It is easy to see that M* 
coincides with the transfer matrix of a rectangular barrier potential flTTll which admits a simple 
explicit expression [2]- It is convenient to write it in the form 



M^o(a) ;= 


M* = A^(nfc^L,n+,n_), 


(26) 


where n± := (n ± n ^)/2, and for all a, b, c G C, 
A^(a, b, c) : = 


(cos a + z b sin a) e 

ikzL 


—ikzL 


I c Sin a e 


—ikzL 


—% csin ae 


(cos a — ib sin a) e 


xkz L 


(27) 


Furthermore, in view of the composition property of transfer matrices, and Eqs. CB. 
and (1271) . we have the following remarkably simple result. 




(28) 


where n™ := (fi^ ± Using a similar analysis we can compute the transfer matrix for 

the TM waves. This yields 

M™ = A^(hfc^L,n™,n™), (29) 

where n™ := (h^ ±d^)/2eh. Equations (1271) - (|29l) imply det M™ = det = 1. Hence we have 
reciprocal transmission; = T’’ =: T. 


6 










For a Veselago slab, which we denote by iS„, we have e = /t = —1, n^ = n^ = l, and (1^ and 
(EHl) give 


M™ = M™ = M„ := 


^—2ikzL g 

Q ^2ihzL 


(30) 


Therefore, in view of is omnidirectionally reflectionless for waves of arbitrary wavennmber. 

Bnt it is not generally invisible, becanse it indnces a phase shift in the transmitted wave; T = 

Let us also note that according to fl27D . (1251) . and (I2U]) . <S is reflectionless for TE/TM modes 
whenever = g. This corresponds to the Brewster angles Of, give by 


cos Ob 


for TE waves, 
for TM waves. 


(31) 


which has the same form for S, i.e., the Brewster angles for metamaterials coincide with those of 
the usual material. 

Another remarkable property of a Veselago slab is that its transfer matrix and consequently its 
scattering features are invariant under space translations and reflections about any plane parallel 

T 

to the slab. This follows from the fact that under a space translation, 2 ; —^ z — a, and the space 

'P 

reflection (parity), z — > —z, transfer matrices transform as 

M ^ T-^M T„, M ^ criM~Vi, (32) 


where 


_ ^iakzcra 


T ‘ 

a • ^ 


^idkz 

0 


■ 0 

1 ■ 


■ 1 

0 

0 

g idkz 

, CTi = 

_ 1 

0 

, era = 

0 

-1 _ 


(33) 


Next, consider placing an arbitrary inhnite planar slab S between a pair of parallel Veselago slabs 
of thickness L. The space-translation invariance of Veselago slabs implies that the transfer matrix 
of both of the Veselago slabs have the form fl30|) . Therefore, if we use M 5 to denote the transfer 
matrix of S, then according to the composition property of the transfer matrices, the transfer matrix 
of the whole system of three slabs has the form 


M = 


_ Ri^RrjTs) R^s/Ts 
-R^s/Ts _ ■ 


(34) 


Here we have made use of ([T9]), assumed that detM^ = 1 (so that the left and right transmission 
amplitudes of S coincide), and denoted the reflection and transmission amplitudes of S by itg and 
Ts- Equation (13T|) shows that no matter what the distance between the slabs are, the transmission 
amplitude of S undergoes a phase shift, according to T —)> , while its reflection amplitudes 

do not change. This marks the application of the Veselago slabs as omnidirectional phase shifters. 


4 Laser Threshold and CPA Conditions 

The condition that an optical potential supports a spectral singularity [22] corresponds to what is 
known as the laser threshold condition in optics [7]. Because the time-reversal of this condition 
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yields coherent perfect absorption (CPA) or antilasing [8l|9], the knowledge of spectral singnlarities 
allows for the determination of the laser threshold and CPA conditions in a variety of optical setnps 
[231 El]- Reference [23] makes use of this approach to provide a detailed description of the behavior 
of the TE and TM waves interacting with a nonmagnetic optically active slab of ordinary material. 
Because spectral singularities (respectively their time-reversal) correspond to the real zeros of the 
M 22 (respectively Mu) entry of the transfer matrix in the complex fc-space [221191123], we can easily 
use fl271) . fl28ll and fl 2 ^ to do the same for the metamaterial slab S. As a result we hnd the following 
compact expression for the laser threshold and CPA conditions. 


^2ikzLn 


TE/TM , 1 

n+ ' -hi 

TE/TM 7 

n+ ' - 1 


7 

5 


(35) 


where 


J 1 for lasing, 
\ -1 for CPA. 


In order to derive a formula for the threshold gain (loss) coefficient and the phase condition for 
lasing (CPA), we introduce 


Si := Re(d), 82 := Im(£), pi := Re(/i), /X 2 := Im(/i), 


so that 


8 — £1 -|- i 82 i h — hi A fp-Sj 


and conhne our attention to situations where the ratio of \ 82 \ and |p 2 | to l^ipi — 1 | are several orders 
of magnitude smaller than 1 , i.e., 

^ 2 ! + |h2 
kihi -1 


In this case, except for the incidence angles close to the Brewster’s angle, we can ignore the quadratic 
and higher order terms in 82 and p 2 , [25]. For example, in view of the fact that = ifi, we can 
write 


n 


yTi/ii + -(£ih2 + hi^2) 


— ^vTi/ii + - (7ilh2 + Ihik2), 


(36) 


where 


J —1 for £1 > 0 and > 0 , 
\ 1 for £1 < 0 and < 0 , 


and we have taken note of the fact that for a metamaterial with £1 < 0 and pi < 0, Re(n) < 0, [3]. 

Observe that fl36|) is consistent with the fact that for a gain metamaterial, 62 < 0 and /i 2 < 0. 
To see this, consider a normally incident TE wave, where 9 = 0, kz = k, and fi = n. Then, as the 
wave travels a distance d inside the slab, the intensity of the electric held, namely \S’{z)\‘^, changes 
by a factor of The gain/loss coefficient is identihed with the quantity g that allows us to 

write this factor in the form Therefore, by virtue of flHB]) . it has the form 


g := -2/cIm(n) = -/c(|£i|/i 2 + Ihik 2 ) for eipi > 0. 


( 37 ) 


Clearly, if 62 < 0 and /i 2 < 0, (/ > 0 and the (meta)material is a gain medium [23]. 








Next, we substitute fl36l) in ffT^ and use ffT6l) and fl37|) to obtain 


fi sec6*(,^ vi + 1 ^ 2 ) 


( 38 ) 


where 



V 2 ■ = 


2 


(kl|h2 + Ihlk2) 


2 k ’ 



(39) 


Inserting (138|) in the expression for and employing the result together with fl3^ and (138|) in 

(153]) . we hnd the following laser threshold (CPA) and phase conditions 
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k 


2'j{nm + 0 ) 
uiL 


(40) 


where we have only retained the leading order terms in powers of £2 and fi 2 and introduced 


rj : = 


0 


J sec 6^ z/i/|/ii| for TE waves, 
\ sec9vi/\ei\ for TM waves. 






I!±l\ 
T]-! J 



for TE waves. 


77^ — 1 


27rm 


In 


f ii±i 

\V-l 



for TM waves. 


(41) 

(42) 


and a mode number m that takes positive integer values. Relations 0391) - 042p show that the sign 
of the real part of the permittivity and permeability does not enter the laser threshold and CPA 
conditions. This shows that as far as these conditions are concerned active metamaterials behave 
exactly like the usual active material. 


5 Metamaterials from Point Interactions and 
Phase-Conjugation 


Consider the matching matrix fl23D that represents the boundary condition at the interfaces of the 
metamaterial slab S and vacuum. We can obtain the matching matrix for the slab S by 

replacing by in fl23l) where 

d™ = 1/d™ = h = l/Uil + i/i2, d™ = 1/d™ = e = |£i| + ie2. (43) 


It is not difficult to see that 

where 


gTE/TM _ gTE/TMgTE/TM 


gTE/TM 


^0 


1 0 
0 — 


(44) 


(45) 


and ^pzQ is the principal argument (phase angle) of i.e. = q;™^™/|q;™'^™|. In other 


words. 


9^0 —^L — 


arctan(p 2 /hi) for TE waves, 
arctan(e 2 /£i) for TM waves. 
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Equation suggests that if we view and as defining certain point interactions, 

then we can identify the transfer matrix for the point interaction given by with the 

product of those associated with and Denoting these with and 

we have 

jyjTE/TM ^ ]yjTE/TM ^TE/TM_ 

In light of ([22]), has the following explicit expression. 


_^TE/TM 


= 


—i sin If e ^Qg 

e^^kzzo (.Qg ip gjj^ p, 


For a Veselago slab, e = /I = 1, (p = 0, 

jyjTE/TM ^ _y^TE/TM 


0 


^-2ikzZo 


^2ikzZo 


0 




(46) 


(47) 


and the transfer matrix of the slab takes the form Substituting (147|) 

in this equations we recover fl30ll . 

Because the bulk transfer matrix M* for our metamaterial slab S is the same as the one for 
S, this suggests that we can view 5 as 5 supplemented with a pair of point interactions given by 
gTE/TM boundaries. In particular, a Veselago slab of thickness L is equivalent to a pair of 

point interactions with center zq = 0,L and matching conditions g™/™ = cr^. Therefore all the 
surprising properties of the metematerial slabs may be realized using certain surface interactions 
that take place at their boundaries. The same conclusion is reached by the authors of Ref. m 
who argue that the properties of a Veselago slab are common to those of a pair of parallel phase- 
conjugating plates. These can also be modeled using an antilinear point interaction of the form 


’ ^(4) 




(Tg ★ 


Y'{z^) _ ’ 


(48) 


where ★ stands for the antilinear operator of complex-conjugation of matrices; given a 2 x 2 matrix 
A and a column vector (2x1 matrix) v, we have -kAv = A* ★ v = A*v*. 

We can purse the approach of [19] to determine the transfer matrix of flTHl) . This results in an 
antilinear matrix of the form 

zo ■= = T_2 o*, (49) 

where N 2 Q and Tq are respectively given by fl^ and fl33|) . It is easy to show that coincides with 
the identity matrix. This is a manifestation of the fact that an overlapping pair of phase-conjugating 
plates should behave exactly like vacuum. 

The composition property of transfer matrices implies that the transfer matrix of a pair of 
parallel phase-conjugating plates located at z = 0 and z = L, is given by 


^0 — ★ Tq^ — — 


^-2ikzL g 

g ^2ikzL ’ 


(50) 


where we have made use of fl49l) . The fact that the right-hand side of this relation coincides with 
the transfer matrix of Veselago’s slab fISUD provides a straightforward proof of the equivalence of the 
latter with a pair of parallel phase-conjugating plates mi- 
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6 Metamaterials and PT-symmetry 

The point interactions defined by fl45l) correspond to the matching conditions of the form 


1p{z„ ) = l/’W), 




( 51 ) 


(62) 


where 

4 2 

It is qnite remarkable that these matching conditions also appear in the study of a class of VT- 
symmetric spectral problems defined on the wedge-shaped complex contours of the form HI: 


r ;= {3 — ztan'd I 3 I | 3 G M}. (53) 

See also FiglU It turns out that the solution of the time-independent Schrodinger equation, 




defined in L‘^{T) is equivalent to the spectral problem defined in L^(M) by the differential equation 
g 2 *,?sgnP) + v{z)ij{z)] = Eij{z) Z e M \ {0}, (54) 


and the matching conditions (l5TD . where v{z) is determined by 1 /( 3 ) and ■d, [H]. 

Let us consider a volume of infinite extent that is separated into two regions by the z = 0 plane 
and filled with nonmagnetic (meta)materials (/i = po) such that the permittivity of the two regions 
differ only by a phase factor according to 


J Sq E- for z < 0, 
\ for z >0, 


where G M’*' and "d G [0, vr/4]. See Fig. |3l 

It is not difficult to see that the TM waves propagating along the ^(-axis in this volume are given 
by d?]) provided that we set d = 0 and demand that J^{z) satisfies fl5T]) and fl54|) with v{z) = 0 
and E := In view of the results of [14], this identifies the behavior of these waves with 

that of the eigenfunctions of the PT-symmetric Hamiltonian operator for a free particle moving on 
the contour fl5^ . The special case £_ = 1 and "d = 7 r /4 corresponds to a volume half-filled with a 
nonmagnetic metamaterial with permittivity —1. In this case J>^{z) solves the spectral problem; 


with 


-sgn{z)—t/j{z) = X^t/jiz), z^O, 

(56) 

^(0+) = ^/^(O"), ^^'(0+) = -^/^'(0“), 

(57) 

= —iE = 

(58) 


Next, consider confining this system to the region defined by 0 < x < L and —L/2 < z < L/2. 
This corresponds to waves propagating in a rectangular waveguide aligned along the y-axis. See 
Fig. m If we suppose that the guide has perfectly reflecting walls and try to obtain its transverse 
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Figure 3: (Color online) An infinite volume of space. The blue region corresponding to z < 0 consists 
of nonmagnetic ordinary material. The pink region corresponding to 2 ; > 0 contains nonmagnetic 
(meta) material. 


modes, we should solve the spectral problem obtained by supplementing fl571) with the Dirichlet 
boundary conditions a.t z = ±T/ 2 , i.e., set 


V'(±L/2) = 0. (59) 

This determines the ^-dependent part of the corresponding standing waves for the magnetic held 
H. The x-dependent part is determined by the solutions of 

- ^T(a^) = Kx{x), x( 0 ) = x(L) = 0 , (60) 

in [0, L]. In other words, the standing magnetic waves for this system are given by the solutions of 
the spectral problem: 

-sgn{z)V'^f{x,z) = Xf{x,z), a: e (0,L), 2 ; e (-f,0) U (0, f), (61) 

= f{L,z) = 0, 2 e[—|,f], (62) 

f{x,0^) - f{x,0~) = dj{x,0^) + dj{x,0~) = f{x,±^) =0, xe[0,L], (63) 

which is dehned on the square □ := { (x, z) | 0 < x < L, \z\ < }. 


The authors of Ref. [12] consider the purely mathematical problem of whether fl 6 T]) - fl 6 ^ can be 
cast into a spectral problem for a genuine self-adjoint operator acting in L^(n). They answer this 
question in the affirmative and show that 0 is an inhnitely-degenerate eigenvalue of this operator and 
that the other eigenvalues are located symmetrically about 0. They also use the separability of flM]) 
to determine the general form of the separable solutions, i.e., the standing waves of our waveguide 
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Figure 4: Cross section of a rectangular waveguide aligned along the |/-axis and half-filled with a 
nonmagnetic metamaterial with permittivity —(pink region). 


in terms of products of the solutions of (1571) and ([60]). The latter gives x(a;) = sin( 7 rnx/L) with 
n = 1 , 2 , 3, • • •, while the former can be expressed in terms of sin and sinh functions [T2] . 

As we noted above (1571) is a special case of the PT-symmetric spectral problem dehned by fl54l) . 
(HJ. and ' 0 (±-|) = 0 , which is treated in [T3|. As shown in this reference, for > 7 r/ 4 , none of 
the corresponding eigenvalues is real. The PT-symmetry of the problem implies that they appear 
in complex-conjugate pairs. For ■d = 7 r /4 which is relevant to our model, Pq = 0 is the only real 
eigenvalue and again the rest of the eigenvalues form complex-conjugate pairs. By analogy to its 
one-dimensional counterpart, we can also reveal the PT-symmetric nature of fl 6 T|) - (|63|). We can 
certainly write the differential equation give in fl 6 T|) as Lf = Xf where 

L := —i sgn{z)V‘^, A = iX. (64) 

Clearly L is PT-symmetric, for Vf{x,z) := f{x,—z) and Tf{x,z) := f{x,z)*. Furthermore, PT 
leaves the set of boundary conditions given in fl62|) and fl 6 ^ invariant. Therefore, written in terms 
of L and A, fl 6 T]) - fl 6 ^ is essentially a PT-symmetric spectral problem. The fact that the spectrum 
of fl 6 T|) - fl63|) is real and symmetric about zero is equivalent to the condition that the spectrum 
of Lf = Xf dehned via the same boundary conditions is purely imaginary and symmetric about 
the real axis. In other words the complex eigenvalues come in complex-conjugate pairs. This a 
characteristic feature of all PT-symmetric spectral problems. 

This example also reminds us of the following rather obvious but mostly unappreciated conse¬ 
quence of the characterization theorems of Refs. [TT] [28] . 

Theorem: Let A he a Hermitian (self-adjoint) operator acting in a separable Hilbert space and 
having a discrete spectrum. Then the anti-Hermitian operator A := iA is pseudo-Hermitian 
provided that the spectrum of A is symmetric about 0 and the eigenvalues of A that are related 
by the reflection about 0 have the same multiplicity. 
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By virtue of a result of [2H] , we can further infer that the operator A commutes with an antilinear 
involution, i.e., it possesses a generalized PT-symmetry [29]. Because the spectrum is not real, this 
symmetry is always broken. 

7 Concluding Remarks 

During the past hfteen years or so, the purely theoretical results of Veselago on unusual properties 
of metamaterials with negative permittivity and permeability have found experimental verihcations. 
The efforts in this direction together with the seminal works of Pendry and his collaborators 
on the application of negative-index metamaterials in constructing perfect lenses and invisibility 
cloaks have made their study a major area of research in theoretical and experimental physics. 

In this article we have employed the idea of point interactions together with the transfer matrix 
formulation of scattering theory to explore some of the basic properties of metamaterials. In partic¬ 
ular, we have computed the transfer matrix of a general (meta)material slab both for TE and TM 
waves and used the result to demonstrate the omnidirectional reflectionlessness of the Veselago slab 
for waves of arbitrary wavenumber. We have shown that sandwiching an arbitrary slab between 
a pairs of Veselago slabs does not affect the scattering features of the original slab except that it 
shifts the phase of the transmitted waves. This marks the application of Veselago slabs as phase 
shifters. We have also established the fact that the Brewster’s angles and the laser threshold and 
CPA conditions do not depend on the sign of the real part of the permittivity and permeability of 
a (meta)material slab, if these have the same sign. 

Another result we have reported in this article is a simple proof of the equivalence of a Veselago 
slab with a pair of parallel phase-conjugating plates m- This required the introduction of a class 
of antilinear point interactions. The transfer matrix of these interactions turn out to be antilinear 
matrices squaring to the identity matrix. These curious observations call for a comprehensive study 
of antilinear point interactions and their physical applications. 

Finally, we have shown how certain conhgurations of metamaterials are related to a class of 
exotic PT-symmetric models dehned on wedge-shaped complex contours. These were originally 
introduced for purely theoretical reasons [H]. Our results provide a simple physical interpretation 
for these models. 
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active metamaterials and helping me hnd and correct typos in the hrst draft of this article. This work 
has been supported by the Scientihc and Technological Research Council of Turkey (TUBITAK) in 
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